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The effect of clouds on the dynamical and chemical evolution of gas-rich 
dwarf galaxies 
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We study the effects of clouds on the dynamical and chemical evolution of gas-rich dwarf galaxies, in particular focusing 
on two model galaxies similar to IZwl8 and NGC1569. We consider both scenarios, clouds put at the beginning of 
the simulation and continuously created infalling ones. Due to dynamical processes and thermal evaporation, the clouds 
survive only a few tens of Myr, but during this time they act as an additional cooling agent and the internal energy of cloudy 
models is typically reduced by 20 - 40% in comparison with models without clouds. The clouds delay the development 
of large-scale outflows, therefore helping to retain a larger amount of gas inside the galaxy. However, especially in models 
with continuous creation of infalling clouds, their bullet effect can pierce the expanding supershell and create holes through 
which the superbubble can vent freshly produced metals. Moreover, assuming a pristine chemical composition for the 
clouds, their interaction with the superbubble dilutes the gas, reducing the metallicity (by up to ~ 0.4 dex) with respect to 
the one attained by diffuse models. 
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1 Introduction 

Intense episodes of star formation (SF) in dwarf galaxies 
(DGs) are associated to the development of large-scale out- 
flows. These have enormous multiple effects on the evolu- 
tion of DGs: the hot gas carries off a large amount of freshly 
produced metals, thermalizes the interstellar medium (ISM), 
reduces the gas content, ceases the SF, and further more. 
The study of this phenomenon through numerical simula- 
tions has therefore been in the focus of several authors in the 
past. The overall picture is that the occurrence of outflows is 
initially driven by the thermal pressure of a hot, high pres- 
surized gas and is favored by a flat distribution of the ISM, 
which allows an easy vertical transport of material, whereas 
horizontal transport is very limited. Therefore, outflows are 
able to eject a large fraction of freshly produced metals but 
only a very limited fraction of ISM (D'Ercole & Brighenti 
1999; Recchi, Matteucci & D'Ercole 2001). 

Most of these studies, however, consider a single-phase 
ISM, although attempts to perform multiphase hydrodynam- 
ical simulations have been made in the past, particularly us- 
ing the so-called chemodynamical approach (Samland, Hensler 
& Theis 1997; Hensler 2003; Hensler, Theis & Gallagher 
2004). In two previous papers, we have studied the dynam- 
ical and chemical evolution of model galaxies similar to 
IZwl8 (Recchi et al. 2004, hereafter Paper I) and NGC1569 
(Recchi et al. 2006). In this work we simulate the same two 
galaxies, but we increase arbitrarily the gas density of some 
specific regions of the computational grid, in order to create 
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a "cloudy" phase, and we address the question how and to 
which extent this new component alters the former results. 

It is particularly interesting to compare the chemical evo- 
lution of models with and without clouds. In fact, clouds 
are expected to dilute the hot metal-rich gas through evap- 
oration, allowing for a reduction of the metallicity without 
altering the abundance ratios (Koppen & Hensler 2005). On 
the other hand, the clouds can act as a cap and hamper the 
development of a galactic wind by means of their drag and, 
if they are dispersed, due to Kelvin-Helmholtz instability 
and evaporation of mass load to the hot outflow. As a re- 
sult these should lead to an increase of the metallicity of the 
galactic ISM, since, as we have seen, galactic winds oth- 
erwise carry metal-enriched gas away from the galaxy. In 
this case, the abundance ratios are affected if the ejection 
efficiencies depend on the different chemical species. 

2 The model 

The simulations are performed by means of a 2-D chemody- 
namical code described in Paper I and references therein. It 
assumes cylindrical symmetry and follows in detail the re- 
lease of mass, individual elements and energy from SNela, 
SNell and stellar winds. As mentioned in the introduction, 
the initial set-up of the models is aimed at reproducing the 
main characteristics of two well-studied gas-rich dwarf galax- 
ies, namely IZwl8 and NGC1569. 

The cloudy phase is introduced by randomly choosing 
grid cells and increasing artificially the density therein and 
in the neighboring cells in order to reproduce a cloud den- 
sity profile similar to the observed ones (namely p c \ oc R~, , 
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Fig. 1 Density contours of the gas for the static cloudy 
model at 4 evolutionary times (labelled, in Myr, on the top 
of each panel). The logarithmic density scale (in g cm~ 3 ) is 
on the right-hand strip. For reference, also the density con- 
tours of the correspondent diffuse model at t ~ 120 Myr is 
displayed (rightmost panel). 
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Fig. 2 Density contours for the infall cloudy model at 9 
evolutionary times (labeled, in Myr, at the top right corner 
of each panel). The logarithmic density scale (in g cm~ 3 ) is 
on the right-hand strip. 



where i? c j is the distance from the center of the cloud; de 
Heij, Braun & Burton 2002). We assume that the chemical 
composition of the clouds is pristine (i.e. without metals). 
The clouds can be either added at the beginning of the sim- 
ulation (static cloudy models) or continuously created dur- 
ing the evolution of the model. In the latter case, we will 
give also an infall velocity of 10 km s _1 to the clouds [in- 
fall cloudy models). We will also consider as comparison 
models obtained without a cloudy phase (diffuse models). 

Since saturated heat conduction is included (Cowie & 
McKee 1977) in the static cloudy model the clouds can only 
act until they vanish due to evaporation, while in the infall 
cloudy model clouds are reborn on timescales of a few Myr. 

3 Results 

3.1 Dynamical evolution 

The snapshots of the evolution of a prototypical static cloudy 
model are shown in Fig. Q] (first 4 panels) together with the 
density distribution on the analogous diffuse model at t ~ 
120 Myr (rightmost panel). This model is aimed at repro- 
ducing IZwl8 and has a SF history characterized by a mild, 
long-lasting episode of SF and a recent burst (see Paper I; 
Recchi & Hensler 2007 for details). The clouds are soon en- 
gulfed in the hot cavity and evaporate on a short timescale 
of a few tens of Myrs. The presence of the clouds is how- 
ever sufficient to create large shears and eddies. Due to these 
structures and due to the evaporation of the clouds within 
the superbubble, radiative losses in this model are signifi- 
cantly larger than in the equivalent diffuse model (by ~ 20 
%). 



At the beginning therefore the evolution of the super- 
bubble is slightly slower than in diffuse models. In particu- 
lar, a galactic outflow starts ~ 75 Myr after the beginning of 
the SF; ~ 5-10 Myr later than the analogous diffuse model. 
However, the presence of clouds strongly distorts the shape 
of the supershell and the highly pressurized gas inside the 
cavity can more easily find regions of lower pressure, from 
which it is easy to pierce the shell and break out. This cre- 
ates the tongues visible in Fig.[T] From these tongues freshly 
produced metals can leave the galaxy. Therefore, in spite of 
the (slightly) slower development of the galactic wind, the 
static cloudy model maintains a large ejection efficiency of 
newly synthesized metals. 

A similar behavior is seen in the infall cloudy model 
(Fig.|2]i. This model has the same setup as the static cloudy 
model seen above, the only difference being the treatment 
of clouds. Bow shocks created by the clouds while falling 
towards the disk are evident in this figure. The interaction 
between clouds and ISM produces also Kelvin-Helmholtz 
instabilities and evaporation, therefore, also in this model, 
clouds do not survive more than a few tens of Myr. The ra- 
diative losses are more significant than in the static cloudy 
model because continuously formed clouds deliver more 
mass to the outflow. At the end of the simulation the ther- 
mal budget of the system is reduced by ~ 35% compared to 
the analogous diffuse model. Consequently, the first weak 
signs of a galactic outflow appear only ~ 90 Myr after the 
beginning of the SF; ~ 20-25 Myr later than in the diffuse 
model. However, analogously to the case of the static cloudy 
model, the ram pressure of the clouds pierces the supershell 
and creates holes through which hot cavity gas (and freshly 
produced metals) can leak out. 
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3.2 Chemical evolution 

The evolution of the mean oxygen abundance in the ISM 
with time for the static cloudy model and the infall cloudy 
model is shown in Fig. [3] in comparison with the evolu- 
tion of the reference diffuse model. The interaction of the 
clouds with the expanding supershell has a non-steady be- 
havior, resulting in moments in which the attained average 
metallicity of the galaxy is larger than in the diffuse model 
and moments in which it is lower. As we have explained 
in the Introduction, this depends essentially on the compe- 
tition between the dilution effect (due to the evaporation of 
the metal-free clouds, which tend to reduce the metallicity) 
and the drag effect (the delayed production of galactic out- 
flows, which tends to increase the metallicity). 

At the end of the simulation, however, the metallicity of 
the cloudy models is lower than the one of the diffuse model 
(by ~ 0.1 dex for the static cloudy model; ~ 0.4 dex for 
the infall cloudy one). This shows that over the long term, 
the dilution effect prevails over the restraint effect of the 
metal-enriched gas and the clouds have on average the ef- 
fect of reducing the metallicity. This is caused by two facts: 
on the one hand, (as explained in Sect. 13. U the interaction 
clouds-supershell distorts the supershell and produces fin- 
gers and tongues through which the metals present in the 
cavity can leave the galaxy, therefore the ejection efficiency 
of metals is not significantly reduced (in some case it is even 
larger) compared to the diffuse model. The second effect of 
metal dilution due to cloud evaporation becomes discernible 
by the tendency that in the infall cloudy model the oxygen 
abundance is continuously reduced with time after a maxi- 
mum. 

Models reproducing NGC1569 show the same kind of 
behavior explained above, therefore we do not describe them 
in this contribution. A comprehensive description of these 
models (and of others in which the IMF and the nucleosyn- 
thetic yields are varied) can be found in Recchi et al. (2006). 

4 Conclusions 

Although models of DG winds show strikingly that the low 
effective yield causes the low metallicity in dlrrs (Garnett 
2002), nevertheless, it is obvious that a vertically extended 
gas disk and/or an enveloping gas reservoir is able to ham- 
per the outflow of hot metal-rich supernova gas. The pre- 
sented results demonstrate that the effects of existing em- 
bedded gas clouds and their evaporation on the effective 
yield is not straightforward predictable but depends on the 
porosity of the gas envelop and its reservoir, i.e. on the de- 
livery of clouds, and must be studied in detail. 

Due to dynamical processes and thermal evaporation, 
the clouds survive only a few tens of Myr. Due to the ad- 
ditional cooling agent, the energy of cloudy models is re- 
duced by 20-40% compared with diffuse ones. The clouds 
delay the development of large-scale outflows, helping in 
retaining a larger amount of gas inside the galaxy. However, 
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Fig. 3 Evolution of 12 + log (O/H) for the static cloudy 
model (dotted line), the infall cloudy model (dashed line) 
and the reference diffuse model (solid line). 

their bullet effect pierces the expanding supershell and cre- 
ates holes through which the superbubble can vent out the 
freshly produced metals. The resulting final metallicity is 
generally smaller (by ~ 0.4 dex at most, depending on the 
model) than the one attained by diffuse models. 

Nevertheless, models with self-regulated SF and self- 
consistent treatment of the ISM, in particular with respect to 
cloudy and hot phases, are requested to evaluate the chemi- 
cal abundances reliably. 
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